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ABSTRACT 

The wind of an asymptotic-giant-branch stars is sufficiently strong that if it is shghtly 
asymmetric, it can propel the star outside of the open cluster of its birth or significantly 
alter its trajectory through a globular cluster; therefore, if these stellar winds are 
asymmetric, one would expect a deficit of white dwarfs of all ages in open clusters and 
for young white dwarfs to be less radially concentrated than either their progenitors or 
older white dwarfs in globular clusters. This latter effect has recently been observed. 
Hence, detailed studies of the radial distribution of young white dwarfs in globular 
clusters could provide a unique probe of mass loss on the asymptotic giant branch 
and during the formation of planetary nebulae both as a function of metallicity and a 
limited range of stellar mass. 
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1 INTRODUCTION 



[Spruit ( 1998[ ) proposed that white dwarfs can acquire their 
observed rotation rates from mild kicks generated by asym- 
metric winds toward the end of their time on the asymptotic 



giant branch (AGB) ([ Vassiliadis fc Wood|p93 l 
|et al 



Fellhauer 



( 2003 1 envoked these mild kicks to explain a putative 



dearth of white dwarfs in open clusters (e.g. Weidemann 



1977 Kalirai et al. 20011. Unfortunately from the dynam- 



ics of open clusters it is difficult to probe the kicks because 
the white dwarfs simply leave the cluster. Although one can 
probe the dynamics of AGB winds directly through observa- 
tions of masers, attempts to look for asymmetries are dogged 
by the variability of the star itself which makes it difficult 
to constrain any reiative motion between the centre of mass 
of the wind and that of the star. The dynamics of globular 
clusters provide a unique environment to probe white dwarf 
kicks because both the magnitude of the expected kick and 
the velocity dispersion of the giants within the cluster are 
on the order of kilometers per second, so one would expect a 
significant signal. On the other hand the escape velocity of 
the cluster may be several times larger so most of the white 
dwarfs remain in the cluster to measure after the kick. The 
expected signature of white dwarf kicks has been observed 
in M4 and NGC 6397 ( [Davis et al.p006 l. 



This letter will examine how white dwarf kicks affect 
the radial distribution of young white dwarfs in globular 
clusters through analytic and Monte Carlo calculations of 
the phase-space distribution function of stars in the cluster. 



2 CALCULATIONS 

Clusters of stars can typically be modelled with a lowered 
isothermal profile (or King model) ( Michie|1963 King|1966 1 
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where e — ~ t^v , is, the gravitational potential, tr is a 
characteristic velocity dispersion and pi is a characteristic 
density. Integrating the distribution function over velocity 
yields the density distribution 
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T he density is typically constant within the core radius, ro = 
•^90-2/(471 po) where po is the central density of the cluster. 
The gravitational potential can be solved self-consistently 
with Eq. Q ( [Binney fc Tremaine|1987| ) using 
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to give a model for the cluster. Because the distribution func- 
tion depends only on constants of the motion (the energy), 
it is constant in time as well. 

With time the kinetic energy within the cluster ap- 
proaches equipartition between the various stars such that 



miCTj = mj(jj (Spitzer 19871. The progenitors of young 



white dwarfs will be the most massive main-sequence stars 
in a cluster at the time, so they will typically have ctto < fc, 
where ctc is the mean velocity dispersion of the cluster; fur- 
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thermore, because these progenitors will only have a small 
fraction of the mass of the cluster, they can be considered as 
massless tracers; their phase space density will be given by 
Eq. ^ with a = aro- The gravitational potential and core 
radius will be determined by Eq. Q-Eq. (|3| with a = ac- 
During their time on the AGB, the stars may lose a large 



fraction of their mass suddenly and asymmetrically ( Vassil 



iadis fc WoodlpM] |Spruit|p98} [Fellhauer et"ar]|2003| ) and 
receive an impulsive kick. To model this kick, the original 
phase-space density can be convolved with a kick distributed 
as a Gaussian: 
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Although this convolved distribution function can be ex- 
pressed in closed form (see the appendix), two key results 
are important. First, the final distribution function depends 
on the value of the gravitational potential and the velocity 
separately and not on the conserved energy alone. Second, in 
the limit of large ^ or small ak the distribution approaches 
a lowered isothermal profile with af ~ a^o + (^k ■ 

Because the distribution no longer depends on constants 
of the motion alone, the distribution function itself will de- 
pend on time, so a Monte Carlo realization of the kicked 
distribution function is helpful to make further progress. 
Specifically, for a variety of values of *I'(0), (tto and at, ten 
thousand stars are drawn from Eq. ([l]) and given a three- 
dimensional velocity kick drawn from a Gaussian of width of 
ak- Finally, each star is evolved forward along its orbit to a 
random phase; this yields an estimate of the radial distribu- 
tion of the young white dwarfs in the cluster after the kick as 
well as the fraction of young white dwarfs that might escape 
the cluster. The potential of the cluster is fixed to be the 
solution to Eq. ([3| for cr = 1; this assumes that the evolving 
AGB stars and young white dwarfs make a negligible contri- 
bution to the mass of the cluster and that the observations 
of the young white dwarfs occur within a relaxation time of 
their formation. 

The result of the Monte Carlo realization is a list of 
radial positions at the moment of the kick and at a random 
time later for those stars that cannot escape the cluster. 
The density distribution given by Eq. Q provides a natu- 
ral characterization of the positions of the particles that do 
not manage to escape the cluster. The best distribution is 
obtained by maximizing the Kolmogorov-Smirnov probabil- 
ity (e.g. Press et al. |1992 I that the list of radial positions 
is drawn from Eq. (|2| with a particular value of cr/; in this 
way the escape fraction and the best-fit af provide an ap- 
proximation of the distribution of the young white dwarfs 
(see Fig.[l]). 



3 RESULTS 

The numerical results depicted in Fig.[2]confirm the analytic 
expectations. For small kicks the final distribution is well 
characterized by 



2 2 2 



(5) 

the esti- 
mate used by Davis et al. ( 2006 1, provides a firm lower limit 
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Figure 1. The cumulative radial distribution for ctxo = = 
0.7(Tc with l'(O) = 3o-^ and 6(T^. The black curve follows the initial 
distribution, the red curve follows the phase-mixed distribution 
after the kick. The green curves give the best fitting functions 
of the form Eq. ([Sjl: a model with 'I'(O) = Scr^ yields a escape 
fraction of 25% and erf = 0.86(Tc, and the model with 'I'(O) = 6cr^ 
gives an escape fraction of 4.4% and af = 0.76(Tc. 



on the kick velocity required to puff up a stellar distribu- 
tion within a cluster. A kick of a given size is less effective 
at changing the radial distribution of stars within clusters 
with deeper potential wells than in shallow wells. 

The second diagnostic is the fraction of the stars that es- 
cape the cluster after receiving a kick. The escape fraction is 
approximately proportional to a^/a^Q, the relative change 
in the kinetic energy of the stars due to the kick. Further- 
more, the fraction decreases exponentially with increasing 
values of '5(0) and increases exponetially with increasing 
values of the initial velocity dispersion. Lighter stars are 
much more likely to escape from clusters with more shallow 
potential wells than heavier stars that are initially more cen- 
trally concentrated. More importantly and in contrast with 
the results for open clusters ( Fellhauer et al.||2003 l, the es- 
cape fraction is typically quite small; it did not exceed 25% 
for any of the models considered (not surprisingly the largest 
number of escapees came from a model with $(0) = 4,a^ and 

O-fc = CTTO = O.SCTc). 

Depending on the distribution of stellar masses in the 
cluster the ratio of ctto to tJc may vary. For simplicity, the 
mass function can be assumed to be a power law. 
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and zero elsewhere. This yields a mean mass of 
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Figure 2. The final values of the best fitting cr for the phase- 
mixed stellar distributions as a function of the initial value of 
""To/^'ci ^(0)/'''c ^^'^ <^k/'^TO- From top to bottom the sets of 
curves are for (Tto/o'c = 0.8, 0.7, 0.6 and 0.5. Within each set 
the curves give \I'(0)/o'^ =4,6 and 8 in black, red and blue re- 
spectively. 



If Minin ~ O.IMto and a = —2.25 (a Salpeter initial mass 
function), M ~ 0.23AfTO and aro ~ 0.5(Tc. 

Typically when the line-of-sight velocity dispersion of a 
cluster is determined, the result is dominated by the bright- 
est stars, that is the giants and the stars near the turnoff; 
therefore, the line-of-sight velocity dispersion measured near 
the centre of the cluster is a direct determination of the value 
of (Jto- Typically, crios ctto as ^'(0) cxd. Specifically, 



• (TtO 



exp 



(*(0)/ 



1.45 



(8) 



to within two percent for ^(O) > 2g%q. For 5'(0) = 4(Txo, 
o"ios = 0.89(T. By looking at the distribution of stellar masses 
in the cluster, the value of <Jc can be estimated from o"to as 
in Eq. ([6| and ([7|. As Tab.[l]shows, the line-of-sight velocity 
dispersion of a globular cluster is typically a few kilometers 
per second, similar to the proposed kick velocities; the two 



clusters observed by Davis et al. ( 2006 1 have velocity disper- 
sions of about 4 kms~ , possibly smaller than the expected 
kicks. Furthermore, the effect could be observable in all of 
the clusters listed with the possible exceptions of 47 Tuc, 
Omega Gen and M22. 
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Figure 3. The fraction of stars that escape the cluster after re- 
ceiving a kick. The escape fraction is typically proportional to 
the square of the relative size of the kick, so this dependence has 
been removed. The black, red and blue points give the results for 
\E'(0)/(T^ =4,6 and 8, respectively, for various values of rang- 
ing from O.IcttO to crxOi a^nd values of cr-po from 0.5 to 0.8crc. 
Where the points are absent the escape fraction was less than one 
part in lO**. 



4 CONCLUSIONS 

Studies of the radial distribution of stars near the turnoff 
and young white dwarfs within globular clusters can pro- 
vide a unique handle on the asymmetry of the rapid mass 
loss while stars are on the asymptotic giant branch. In the 



case of an open cluster as studied by Fellhauer et al. (20031, 



the assumed dearth of young white dwarfs could only pro- 
vide an lower limit on the kick; to get a better result requires 
searching for the white dwarfs that have just left the clus- 
ter. In a globular cluster one can get a quantitative estimate 
of the typical kick velocity from the observed distributions 
because very few stars typically escape the cluster; further- 
more, by comparing the distributions in globular clusters of 
different ages and metallicities one can probe the mass loss 
for a variety of stellar masses and chemical compositions. 
By comparing the distribution of white dwarfs of different 
ages up to and greater than the relaxation time for the clus- 
ter, the process of dynamical relaxation of the white dwarf 
population can be constrained observationally. 

Such a study requires observations of many stars 
throughout a globular cluster, specifically the accurate ra- 
dial positions of the brighest main-sequence stars (near the 
turnoff) and the brightest white dwarfs. Because informa- 
tion from the central regions of the cluster is most helpful, 
understanding or at least minimizing the effects of confusion 



is crucial; the space-based observations presented in (Davis 
1,2006) indeed found evidence of white dwarf kicks in 



let al. 



NGC 6397 and M4. Further studies of resolved stellar pop- 
ulations in both open clusters and globular clusters could 
provide further probes of this violent stage in the evolution 
of stars like our Sun. 
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Table 1. Kinetic parameters for a few nearby globular clusters whose radial stellar distribution could be probed with JWST (Richer 
priv. comm.). The concentration or the ratio of the core radi us to the tidal radius (e = login Ttjr^ are from |Harris| | |1996l l , a dash denotes 
that the core has collapsed. The values of >I/(0)/(t^ are from [McLaughlin fc van der Marel| |2005| with the exception of NG C 6837 where 
it was inferred from the concentration. The value of Oc was obtained by assuming that M /M'yq = 0.3/0.8 = 0.375. Refs: (a) |McLaughlin| 
|et al.| | |2006l l; (b) |Reijns et al.| | |2006| ; (c) [Peterson fc Latharn] l |1986| ; (d) |Pryor fc Meylan| (l993| . 



Cluster 


o-los [km/s] 


(Tc [km/s] 


c 


*(0)/<x2 


NGC 104 (47 Tuc) 


11.6±0.8(=') 


19 


2.03 


8.60±0.10 


NGC 5139 (Omega Cen) 


15(b) 


25 


1.61 


6.20±0.20 


NGC 6121 (M4) 


3.9±0.7('=) 


6.4 


1.59 


7.40±0.10 


NGC 6397 


3.5±0.2('^) 


5.7 






NGC 6656 (M22) 


8.5±1.9W 


14 


1.31 


6.50±0.20 


NGC 6752 


4.5±0.5W 


7.3 






NGC 6809 (M55) 


4.2±0.5W 


6.8 


0.76 


4.50±0.10 


NGC 6838 (M71) 


2.8±0.6('=) 


4.6 


1.15 


5.4 



APPENDIX : KICKED PHASE-SPACE 
DISTRIBUTION FUNCTION 

The convolution of the lowered isothermal profile distribu- 
tion function with a Gaussian kick in velocity is given by 



/fin 



2\3/2 



d V exp 



(v- 



2-i 



/initial {x,v'') .(9) 



Because both the kick and basic distribution function are 
Gaussian it is not surprising that the final distribution func- 
tion can be expressed in closed form. If there was no restric- 
tion on the velocities in the distribution function in Eq. ([T]) , 
the convolution would also be a Gaussian with ctj = a^+a^; 
however, the restriction on the velocities complicates the re- 
sult: 



+ 



1 2-Ka exp 

where 

N = 



[g{v)+9{-v)]} (10) 
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(12) 



In the limit of large 'i' or small ak the distribution ap- 
proaches a lowered isothermal profile with a ^ cjf. 
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